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1. Introduction 

A number of recent experimental developments have refocused attention on a relatively 
simple picture [H |2] of the origin of superconductivity in quasi-two-dimensional 
correlated electron compounds. Exchange of spin density wave (SDW) fluctuations 
induces an attractive interaction between electrons in a spin-singlet, even-parity channel, 
with the pair wavefunction changing signs between regions of the Fermi surface 
connected by the SDW ordering wavevector. 

In the pnictide superconductors, experiments identified the nature of the SDW 
order and the configuration of the Fermi surface, and this theory then implies 
superconductivity with s± pairing [21 H], for which strong evidence has recently appeared 

In the cuprates, this picture correctly predicted the d-wave pairing signature. 
However, it faced the difficulty that no strong SDW fluctuations were experimentally 
observed in the optimal hole-doping region, where the superconductivity was the 
strongest. A potential resolution is offered by the theory of competing orders [HI [TJ [8]. 
In the strong-coupling regime, the superconducting pairing amplitude has a significant 
feedback effect on the SDW fiuctuations, and shrinks the region of long-range SDW 
order [9]; consequently the SDW fiuctuations are also substantially reduced at low 
temperatures in the region of strongest superconductivity. An immediate consequence 
of this theory is that the SDW order should re-emerge when superconductivity is 
suppressed by an applied magnetic field. Such a re-emergence has been observed 
in a variety of experiments on both the LSCO and YBCO series of compounds 
[ini im [I2l [ISl E] . The high-field quantum oscillations [15] are also naturally understood 
in this theory, with the re-emergent long-range SDW, or associated, order breaking up 
the Fermi surface into Fermi pockets. 

Motivated by these developments, we recently presented [16] a detailed low energy 
description of the vicinity of the SDW transition in two-dimensional metals. Here we 
will review our main results on the additional instabilities present near such a transition. 
We will phrase our results for the case of the cuprate Fermi surface configuration, 
although there are natural extensions to the pnictides [9]. As expected [H [T7], we 
find an instability to spin-singlet d-wave pairing. However, the pairing instability is 
enhanced from the familiar BCS logarithmic divergence to a logarithm squared. We 
also find a dominant secondary instability to a certain modulated bond order which is 
locally an Ising-nematic order: this instability is also associated with a susceptibility 
which diverges as a logarithm squared, but with a co-efficient which is smaller than that 
of the pairing instability. 

2. Low energy theory 

We begin with a Hubbard-like model for electrons Qo- on the sites (i) of the square 
lattice {a =ti is the spin index) with dispersion eg as a function of wavevector k. We 
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decouple the repulsive interaction in the spin channel by a Hubbard- Stratonovich field 
0" (a = x,y,z). This bosonic field 0" represents collinear SDW order at wavevector 
Q = (tt, tt). In this manner we obtain the familiar "spin-fermion" model for the vicinity 
of the SDW ordering transition [H [T7] with the Lagrangian 



^=/^4.(|: + -.^)---^E<^?4r° 



+ d'x 



^,{drrr + li^.rr + \{rf + \wff 



(1) 



Here r is imaginary time, r"" are the Pauli matrices, and we have partially integrated 
out high energy electrons to generate a bare 0^ field theory Lagrangian for the SDW 
order parameter. The couplings r and u tune the strength of the SDW fiuctuations: in 
mean-field theory there is an onset of long-range SDW order at r = 0. 

The Lagrangian in ([1]) is not yet in the form of a low energy theory amenable to a 
continuum renormalization group (RG) analysis. This is because of the transfer with the 
large wavevector Q in the spin-fermion term A: this will scatter fermions with momenta 
well away from the Fermi surface, and so involves high energy excitations. 

We focus on the low energy sector by zeroing in on the "hot spots". These are 
special locations on the Fermi surface where £^_,_q = £%'■ for these points, both the 
initial and final electron states in the spin-fermion coupling can be right on the Fermi 
surface, and so the scattering from the SDW fiuctuations is the strongest. For the 
electron dispersion appropriate to the cuprates, there are n = 4 pairs of hot spots, as 
shown in Fig. [H We introduce fermion fields {;ip{„, "^D) ^ = l.-.n, a =1^, for each pair of 
hot spots. Lattice rotations map the pairs of hot spots into each other, acting cyclically 
on the index L Moreover, the two hot spots within each pair are related by a refiection 
across a lattice diagonal. It will be useful to promote each field if) to have A^-fiavors with 
an eye to performing a 1/A^ expansion. (Note that in Ref. [17], the total number of hot 
spots 2nN is denoted as A^.) The fiavor index is suppressed in all the expressions. The 
low energy effective theory is then given by the following Lagrangian density in 2+1 
spacetime dimensions for the boson 0" and the fermions ip 

+ Ar (^ilC^L' + ^2lC'^L') • (2) 

The first line in Eq. ([2]) is the usual 0(3) model for the SDW order parameter, the 
second line is the fermion kinetic energy, and the third line is the interaction between 
the SDW order parameter and the fermions at the hot spots. Here, we have linearized the 
fermion dispersion near the hot spots, and v^ are the corresponding Fermi velocities. It 
is convenient to choose coordinate axes along directions x = 4^(1, 1) and y = 4j(— 1, 1), 
so that 

vt^ = {vx,Vy) , iff^ = {-Vx,Vy), (3) 
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Figure 1. Square lattice Brillouin zone sliowing ttie Fermi surface appropriate to 
the cuprates. The filled circles are the hot spots connected by the SDW wavevector 
Q = (tt, tt). The locations of the continuum fermion fields tpl and ^2 ^^^ indicated. 
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Figure 2. Configuration of the £ — 1 pair of hot spots, with the momenta of the 
fermion fields measured from the common hot spot at fc = 0, indicated by the dark filled 
circle. The Fermi velocities ui^2 of the '01.2 fermions are indicated. The momentum 
components of the ipliP) fermion parallel {p«) and orthogonal {p±) to the Fermi surface 
are indicated. 



The velocity ratio a = Vy/v^ plays an important role in the RG analysis of (|2]). The 
Fermi velocities ([3]) are indicated in Fig. [2l The other Fermi velocities are related by 
rotations, v^ = {Rn/2) 



-^if=\ 
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3. Pseudospin symmetry 

We now note a crucial symmetry of the theory ([2]), which will be important for the 
instabilities discussed in this paper. Besides the microscopic translation, point-group, 
spin-rotation and time-reversal symmetries, the low energy theory possesses a set of four 
emergent SU(2) pseudospin symmetries associated with particle- hole transformations. 
Let us introduce a four-component spinor. 

We will denote the particle-hole indices in the four- component spinor by a,f3. The 
spinor (jlD satisfies the hermiticity condition, 

= ^r (5) 




Then, the fermion part of the Lagrangian (12]) can be rewritten as, 

+ -A0" ■ ('^f r"^^ + ¥^t''¥^ . (6) 

Now the Lagrangian ([6]) and the condition (j5]) are manifestly invariant under, 

SU(2),: ¥,-^Ut^ , r^r (7) 

where the U^ are SU(2) matrices. Note that there are 4 independent SU(2) pseudospin 
symmetries, one for each pair of hot spots. The diagonal subgroup of ([7]) is associated 
with independent conservation of the fermion number at each hot spot pair. 

The symmetry ([7]) is a consequence of linearization of the fermion spectrum near the 
hot spots and is broken by higher order terms in the dispersion. The diagonal subgroup 
noted above is preserved by higher order terms in the dispersion, but is broken by four- 
fermion interactions, which map fermion pairs from opposite hot spots into each other. 
Both symmetry breaking effects are irrelevant [16] in the low energy limit used to derive 

This large quadrupled SU(2) symmetry is thus a generic feature of the vicinity of 
the SDW ordering transition in metals. It should be contrasted with the more familiar 
lattice SU(2) pseudospin symmetry of the Hubbard model p^[T9]: there is only a single 
such SU(2) symmetry, present at half-filling, and only if the dispersion has a particle- 
hole symmetric form i.e. the electron hopping is always between opposite sublattices. 
No such restrictions are placed in our case. We will also find a connection between 
instabilities in the particle-hole and particle-particle channels, but both will apply to 
the generic spin- fermion model in (l2]). 

The pseudospin symmetry ([7]) constrains the form of the fermion Green's function 
to be, 

- (^L^ri') = ^''^^^AP^..'G',{^ - x') (8) 
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Figure 3. Pairing of the electrons at the £ = 1,3 hot spots of Fig. [T] Electrons 
at opposite ends of the arrows form spin-singlet pairs. The pairing amplitude is in 
the spin-singlet channel and has opposite signs on the two arrows. There is a similar 
pairing instability between the £ — 2,A hot spots. 



which imphes, 

Gf (x - x') = -G'iix' - x) (9) 

The corresponding expression in momentum space, Gf{k) = —Gl{—k), imphes that 
the location of hot spots in the BriUouin zone is not renormahzed by the spin wave 
fluctuations in the low energy theory. 



4. Pairing instability 

A detailed RG analysis of the SDW transition described by ([2]) was presented in Ref. [16] 
using the framework of the 1/A^ expansion. We will not describe this here, but will use 
a few key results below. 

Let us directly proceed to discuss the nature of the pairing instability. The pairing 
is induced by 0-fluct nations, which play the role of the phonon mode in conventional 
BCS theory. The pairing involves 2 pairs of hot spots, consisting of time-reversed pairs 
of electrons, as shown in Fig. |3l The spin singlet superconducting order parameter, 
associated with Fig. |3] is 

..'(^11^11-41^^11) (10) 



c 



and the relative sign between the two terms is the signature of d-wave pairing. Thus, 
the operator (fTOj) is odd under a reflection about the lattice diagonal P : {k^, ky) — >■ 



^ fCj;, Kyj 
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Figure 4. Cooperon diagram associated with the pairing instability of the £ = 1,3 
hot spots. The full lines are the fermion propagators, as indicated, and the dashed line 
is the propagator of the SDW fluctuation 0". The circled signs indicate the signs of 
the pairing amplitude of the corresponding pairs of fcrmions. 
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Figure 5. The first two terms in the ladder series of Fig.HJ with frequency /momentum 
labels, as evaluated in Eqs. (fT3| and (fT4|) . 



The pairing instability in the channel (ITOl) is described by computing the four- 
fermion scattering amplitude, 

rBCs(fc, wi, -k, W3; k', u[, -k', W3) = 

'^lik',u;[)^U-k',uj',) - ^l{Pk\uj[)ijl,{-Pk\uj',) 



^cra'^pp' 



,lt 



3t 



,lt, 



X [^tl'i-k,u:s)^ll{k,io,)-r2l'{-Pk,uJs)^2l{Pk,io,) 



:ii) 



where it is understood that the external fermion Green's functions have been truncated. 
As in conventional BCS theory, we will study the contribution of the Cooperon ladder 
diagrams illustrated in Fig. HJto the amplitude f lTTl) . Here we will only examine the first 
two terms in the series, as shown in Fig. [51 The dashed line denotes the propagator 
of the field, representing the SDW fluctuations. For this we use the overdamped 
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paramagnon form appropriate to the 1/A^ expansion, 

D{u,q) = ^ I (12) 

where 7 = n\^ / {2TiVxVy) is the damping constant obtained from a decay of a into a 
particle-hole pair. 

The strongest infra-red enhancement of the scattering amplitude by SDW 
fluctuations occurs in the kinematic regime when the incoming fermions have momenta 
at hot spots fc = 0, while the outgoing fermions are "cold", i.e. their momentum k' is 
along the Fermi surface and their frequency satisfies '^\uj\ <^ fc'^. For simplicity, we set 
all the external frequencies to be equal ui = uj^ = uj[ = uj'^ = u . With the above choice 
of external momenta the diagram in Fig. [5^ evaluates to 

19X2 

r*B^- = ^^ (13) 

Next, we consider the diagram for the BCS amplitude in Fig. |5b, which evaluates 
to 

Tbc? = 36A^ j ^k^D{i^^ ^)G\{^r + 00, i)Gl{-ir + 00, -i)D{i,, i- k') 

(14) 

where G denotes the fermion Green's function. The main contribution to Eq. flT^ comes 
from momenta |£| ^ \k'\, 7|^r| ^ ^'^ such that the first scattering in Fig. [5)d is "soft", 
while the second one is "hard" . Then, 

Tl'^^i -^ ^^1 ^^D{U)Gl{(ir + u^A)Gl{-lr + u^.-i) (15) 

N-f\k'^\ J {^T^r 

Now we need the form of G near the SDW critical point. This was described in some 
detail in Ref. [12] to one-loop order in the 1/N expansion: it has a complex structure 
depending upon the value of frequency and momentum from the hot spot. For the BCS 
scattering amplitude, it turns out [16] that the dominant contribution comes from the 
Fermi liquid poles on the "cold lines" in the vicinity of the hot spot. Let us define p± 
as the distance to the Fermi surface, and p\\ as the distance to the hot spot (see Fig. [2]); 
then there are well-defined Landau quasiparticles for p_L ^ py and 7|a;| -C Py, with 

g(^,Pl- . ^^?\ • (16) 

%u -vf[p\\)pi. 

The Fermi velocity vp and the quasiparticle residue Z both vanish linearly with p\\ as 
we approach the hot spot [16], 






where recall a = Vy/vx- We now insert flT2|) . f lT6|) . and flT7|) into flTSj) and obtain 



.Hoop 36A'^ f d£\\ f dir C di\ 1 



''''^ m^\k'^\J 27r J 271 J 27r(7|4|+£j) 
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where £||__l are defined with respect to the 1P2' Fermi surfaces, as in Fig. |2l Recalhng 
that the internal fermions are taken to be cold, i.e. jlirl ^ ^n, we may approximate 
the bosonic propagator in Eq. f lTSjl by its static value. After changing variables to 
e = VF{i\\)£± we then obtain, 

^iioop ^ _ 36A^ r dl\\ Z\i{) 
''''' mj\k'^\J 27rt;^(£|,K| 

r dir f de 1 1 

"" J^ier\<q 27 y 2^ [i{ir + u) - e] [t{ir - u) + e]' 

The integral over ir, e has the form familiar from Fermi-liquid theory and gives the 
usual BCS logarithm, 

dir f de 1 1 _ 1 Afl , . 

27r i 27r[i(4 + w)-e][i(4-w) + e] 27r °^ |u;| ^ ' 

where Ap^ is the frequency/energy cut-off, which in the present case is KpL ~ ^\h- ^^ 
course, for the above form to hold, we need \uj\ <^ ^fl- Thus, 

T.1100P f 12A^ \ 3A2 r'^"' ,, Z\h) , q 

log -1-1 (21) 



where we have cut-off the ultra-violet divergence of the £|| integral by \k'\. Comparing 
Eq. fl2T|) with Eq. (IT3l) we see that there is an enhancement of the Cooperon propagator 
by the factor 

7r(a^ -I- 1) 7|a;| 

which has the promised log-squared form. Note that this is not suppressed by a factor 
of 1/A^. As we noted, the critical SDW fluctuations enhance the BCS logarithm to the 
stronger divergence above. Note that this divergence came from "cold" internal fermion 
lines in Fig. |5]d; thus, similar effects will also be present in higher order graphs in the 
Cooperon ladder. 

It is not clear how to improve Eq. (l22l) using the RG. However, we can note that 
the coupling a is of order unity, and so the pairing is enhanced as the frequency crosses 
the Fermi energy. We also note that in the two-loop RG, the coupling a = Vy/vx has a 
flow towards weak coupling 

da 12 a^ 

m ^ ~TTnN{l + a^) 
but it not appropriate to simply insert the integrated value from this flow into the 
pairing enhancement. 
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Figure 6. Pseudospin partner of the Cooperon susceptibility in Fig |4l obtained by 
reversing the line of the £ = 3 fermion. As described in the text, this represents the 
susceptibility to a modulated bond order which is locally an Ising-nematic order. 



5. Bond order instability 

Let us now apply the pseudospin transformation in ((Tj) to only the i = 3 pair of hot 
spots, while leaving the i = 1 pair unchanged. This will transform the Cooperon 
instability in the particle-particle channel to an instability in the particle-hole channel. 
Diagrammatically, this corresponds to reversing the i = 3 line in Fig. H] to obtain 
the susceptibility in Fig. |6l By the symmetry of the low energy theory in Section |3l 
it would seem we can immediately conclude that the scattering amplitude in Fig. [6] 
should also have a log-squared divergence at one loop order as in (!2T|) . However, as 
was shown in [16], the naively irrelevant curvature of the Fermi surface, which breaks 
the pseudospin symmetry, modifies the result in the particle-hole channel at very low 
frequencies |a;| <^ ^l^^^'P, where A is a co-efficient related to the curvature. It was found 
that in this regime, the log-squared divergence did indeed survive, but with a co-efficient 
which was smaller by a factor of 3. Thus in place of the factor in Eq. ( 122|) for Fig. HI 
the amplitude in Fig. |6] has the enhancement factor 

log^-^. (23) 



1 + 



a 



37r(a2 + l^ 



71(^1 



Moreover, for higher frequencies, u ^ A|A;'|^, the particle-particle and particle-hole 
channels are degenerate. 

What type of ordering does the susceptibility in Fig. E] correspond to? We can 
deduce this by performing the pseudospin transformation on Fig. |3l this changes the 
i = 3 fermions from particles to holes, leading to Fig. [71 Using Eq. f lTU]) . this makes it 
clear that the scattering amplitude in Fig. [6] is that associated with the ordering of the 
operator 



O 



'It 



If, 



^i>L - r2ir2. 



(24) 



and the sign between the two terms is the pseudospin descendant of the d-wave pairing. 
Unlike, the Cooper pair operator, the operator O carries a non-zero momentum Qi, as 
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Figure 7. Spin singlet density operators (~ ip^ijj) of the electrons at the £ = 1,3 hot 
spots of Fig. [1] shown with an arrow pointing from the Brillouin zone location of ip to 
that of i/it. The dashed arrows are the density operators in the first Brillouin zone. 
The full arrows are in an extended zone scheme which shows that these operators have 
net momentum Qi = 2Ky{—l,l), where (K^^Ky) is the location of the i = 1, i = 1 
hot spot (we have K^ + Ky = tt). The density operator with opposite signs on the 
two arrows is enhanced near the SDW critical point. Similarly the £ = 2,4 hot spots 
contribute density operators at Q2 — '^Ky{l, 1). 



shown in Fig. [71 For the cuprate Fermi surface, this is a small momentum oriented 
along a square lattice diagonal. This non-zero momentum would suggest that O is a 
charge density wave at wavevector Qi. However, this does not account for the relative 
sign in i^^, which implies that O changes sign upon reflection about the square lattice 
diagonal parallel to Qi- The situation becomes clearer, when we write (12^ in terms of 
the underlying lattice fermions c^^. The state in which (O) 7^ has 



'"fc-Qi/2,o-'"fcH 



1/2 



(cos kx — cos ky 



(25) 



Despite the (i- wave-like structure, this order is not the popular li-density wave ^20j; the 
latter is odd under time-reversal, and in the present notation takes the form 



'^k-Q/2,a'^k+Q/2,a 



(sin k^ — sin k„ 



(26) 



with Q = (7r,7r). The order in ( l26l) is not enhanced near the SDW critical point, while 
that in (125!) is. By taking the Fourier transform of (!25|l . it is easy to see that O does not 
lead to any modulations in the site charge density (ct^c^o-), and so it is not a charge 

density wave. The non-zero modulations occur in the off-site correlations ( ctg-C^o- ) with 
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Figure 8. Plot of the bond density modulations in ([271) . The lines are the links of the 
underlying square lattice. Each link contains a colored square representing the value 



Csa ) I where f and s are the sites at the ends of the link. We chose the ordering 



of (ci 

wavevector Qi ^ (27r/16)(l, — 1). Notice the local Ising-nematic ordering, and the 
longer wavelength sinusoidal envelope along the diagonal. 



f^ s. For r and s nearest-neighbors, we have 






(0)e 



iQvif+s}/2 



~r C.C. I ^O'l^—s^x ~r Og—j^x ^r~s^y ^s 



'r,y\ 1 



(27) 



where x and y are unit vectors corresponding to the sides of the square lattice unit cell. 
The modulations in the nearest neighbor bond variables (ct^c^+^.o-) and (c^^uCipj^y^fjj 
are plotted in Figs. [S] and O These observables measure spin-singlet correlations across 
a link: if there are 2 electrons on the 2 sites of a link, this observable takes different 
values depending upon whether the electrons are in a spin singlet or a spin triplet state. 
Thus O has the local character of a valence bond solid (VBS) order parameter. The 
first factor on the rhs of Eq. fl27p shows that the VBS order has modulations at the 
wavevector Qi along a square lattice diagonal. As we saw in Fig. [71 |Qi| is small and 
so the first factor in f[27[) contributes a relatively long-wavelength modulation, as is 
evident from Figs. [S and [9l This long- wavelength modulation serves as an envelope to 
the oscillations given by the second factor in ([27[) . The latter indicates that the bond 
order has opposite signs on the x and y directed bonds: this short distance behavior 
corresponds locally to an Ising-nematic order, which is also evident in Figs. [8] and [9] 
The ordering in ( l27l) becomes global Ising-nematic order in the limit Qi -^ 0. It is 
interesting and significant that a Qi = Ising-nematic instability appears in numerical 
weak-coupling renormalization group analyses [2T| [22]. 
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Figure 9. As in Fig. [51 but for orderings afong both Qi — (27r/16)(l, — 1) and 
Q2 = (27r/16)(l,l). 



6. Conclusions 



The ordering found in Section |5] has not (yet) been directly detected in any experiments 
on the cuprates. However, it does have a natural connection to a variety of recent 
experiments [231 1211 [251 ISSl |27]. The ordering is present only in bond observables 
(see Figs. |H]and|ni), and so relates to the bond-centered modulations seen in scanning 
tunnelling microscopy. Also, we have so far performed only a linear stability analysis 
near the SDW critical point. After accounting for non-linearities, and the interplay 
with the superconducting order parameter, generically Qi will move away from the 
value determined by the location of the hot spots: there is no symmetry which pins the 
value of Qi. In general, non-linearities prefer commensurate wavevectors, and so a shift 
to Qi = is not unreasonable. In that case, as we noted above, global Ising nematic 
order would appear. 

Applying a magnetic field will strongly suppress the instability to (i-wave 
superconductivity noted in Section HI However, the particle- hole instability of Section [5] 
occurs in the zero charge channel, and so should be insensitive to the applied field. 
Thus it is possible that the ordering in fl25|) becomes the dominant instability at large 
magnetic fields, and so plays a role in the structure of the Fermi surface and the quantum 
oscillations. 



Instabilities near the onset of spin density wave order in metals 14 

Acknowledgements 

This research was supported by the National Science Foundation under grant DMR- 
0757145, by the FQXi foundation, and by a MURI grant from AFOSR. 

References 

[1] Scalapino D J, Loh E and Hirsch J E 1986 Phys. Rev. B 34, 8190. 

[2] Scalapino D J 2010 A common thread Preprint arXiv:1002.2413 

[3] Mazin I I, Singh D J, Johannes M D and Du M -H 2008 Phys. Rev. Lett. 101, 057003. 

[4] Seo Kangjun, Bernevig B A and Hu Jiangping 2008 Phys. Rev. Lett. 101, 206404. 

[5] Hanaguri T, Niitaka S, Kuroki K and Takagi H 2010 Science 328, 474. 

[6] Sachdev S and Zhang S-C 2002 Science 295, 452. 

[7] Demler E, Sachdev S and Zhang Y 2001 Phys. Rev. Lett. 87, 067202. 

[8] Sachdev S 2010 Physica Status Solidi B 247, 537 {Preprint arXiv:0907.0008) 

[9] Moon E-G and Sachdev S 2010 Quantum critical point shifts under superconductivity: the 
pnictides and the cuprates Preprint arXiv:1005.3312. 
[10] Lake B et al 2002 Nature 415, 299. 
[11] Khaykovich B, Wakimoto S, Birgeneau R J, Kastner M A, Lee Y S, Smeibidl P, Vorderwisch P 

and Yamada K 2005 Phys. Rev. B 71, 220508. 
[12] Chang J et al 2008 Phys. Rev. B 78, 104525. 
[13] Chang J et al 2009 Phys. Rev. Lett. 102, 177006. 
[14] Haug D et al 2009 Phys. Rev. Lett. 103, 017001. 
[15] Doiron-Leyraud N, Proust C, LeBoeuf D, Levallois J, Bonnemaison J-B, Liang R, Bonn D A, 

Hardy W N and Taillefer L 2007 Nature 447, 565. 
[16] Metlitski M and Sachdev S 2010 Phys. Rev. B 82, 075128. 

[17] Abanov Ar, Chubukov A V and Schmalian J 2003 Advances in Physics 52, 119. 
[18] Yang C N 1989 Phys. Rev. Lett. 63, 2144. 
[19] Kee H-Y 2010 Annals of Physics 325, 1260. 

[20] Chakravarty S, Laughlin R B, Morr D K and Nayak C 2001 Phys. Rev. B 63, 094503. 
[21] Halboth C J and Metzner W 2000 Phys. Rev. Lett. 85, 5162. 
[22] Yamase H and Kohno H 2000 J. Phys. Soc. Jpn. 69, 2151. 

[23] Ando Y, Segawa K, Komiya S and Lavrov A N 2002 Phys. Rev. Lett. 88, 137005. 
[24] Kohsaka Y et al 2007 Science 315, 1380. 
[25] Hinkov V, Haug D, Fauque B, Bourges P, Sidis Y, Ivanov A, Bernhard C, Lin C T and Keimer B 

2008 Science 319, 597. 
[26] Daou R et al 2010 Nature 463, 519. 
[27] Lawler M J et al 2010 Nature 466, 347. 



